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We demonstrate numerically that superradiance could play a signifi¬ 
cant role in nonphotochemical quenching (NPQ) in light-harvesting 
complexes. Our model consists of a network of five interconnected 
sites (discrete excitonic states) that are responsible for the NPQ 
mechanism. Damaging and charge transfer states are linked to their 
sinks (independent continuum electron spectra), in which the chemi¬ 
cal reactions occur. The superradiance transition in the charge trans¬ 
fer (or in the damaging) channel, occurs at particular electron trans¬ 
fer rates from the discrete to the continuum electron spectra, and 
can be characterized by a segregation of the imaginary parts of the 
eigenvalues of the effective non-Hermitian Hamiltonian. All five ex¬ 
citonic sites interact with their protein environment that is modeled 
by a random stochastic process. We find the region of parameters in 
which the superradiance transition into the charge transfer channel 
takes place. We demonstrate that this superradiance transition has 
the capability of producing optimal NPQ performance. 

Non-Hermitian Hamiltonian | nonphotochemical quenching | superradiance | 
sink | 

Abbreviations: NQP, nonphotochemical quenching; RTP, random telegraph process; 
LHCs, light-harvesting complexes; ET, electron transfer. 

W hen sunlight intensity is too high, damaging pro¬ 
cesses (such as singlet oxygen production) occur, that 
can destroy the photosynthetic organisms. To survive in¬ 
tense sunlight fluctuations, photosynthetic organisms have 
evolved many protective strategies, including nonphotochem¬ 
ical quenching (NPQ) [T], (See also references therein.) Using 
this strategy, light-harvesting complexes (LHCs) experience 
geometrical reorganizations due to the conformational changes 
of their protein environments. As a result, the damaging chan¬ 
nels are partly suppressed, and the excessive sunlight energy 
is transfered to the non-damaging, quenching channel(s), that 
are opened in this regime. 

In modeling the NPQ mechanisms, it is possible to charac¬ 
terize the sites of the LHCs by the discrete excitonic energy 
states, |n), n being the number of site, associated with the 
light-sensitive chlorophyll or carotenoid molecule. Then, both 
the damaging and undamaging channels can be characterized 
by their corresponding sinks, ISVi), that provide independent 
continuum electron energy spectra [2]. These sinks can have 
very complex structures, and they can be responsible for pri¬ 
mary charge separation processes, as in the photosynthetic 
reaction centers (RCs) [ 3 ] [ 3 ], for creation of charge transfer 
states, for singlet oxygen production [2], an d for many other 
quasi-reversible chemical reactions. Generally, a particular 
sink, | Sri ), is connected to a particular site, |n). The energy 
transfer from this state to its sink is characterized by the elec¬ 
tron transfer (ET) rate, r n . For those sites which are not 
connected to sinks, the corresponding ET rates vanish. Under 
reasonable assumptions, this type of model can be described 
by an effective non-Hermitian Hamiltonian HU SUSIE]. 

Then, this approach becomes similar to those used in de¬ 
scribing the so-called “superradiance transition” (ST) in sys¬ 
tems in which the discrete (intrinsic) energy states interact 
with the continuum spectra 0 0 0 OB [ft] [12]. In these 


systems, the ST usually occurs when the resonances start to 
overlap. Namely, when the half-sum of the imaginary parts of 
the neighboring eigenenergies (decay widths) becomes of the 
order of the spacing between the resonances (distance between 
the real parts of the neighboring eigenenergies). 

For small ET rates, the decay widths are small, and res¬ 
onances are isolated. However, if the r n approach critical 
values, at which the resonances begin to overlap, then a seg¬ 
regation of the widths occurs. In this case, most of the decay 
width is allocated to a few short-lived “superradiant” states, 
while all other states become long-lived (and effectively de¬ 
coupled from the environment). This segregation is called the 
ST. The ST is a quantum coherent effect and, in biological 
systems, it should be considered (see, for example, [ED 01) in 
close relation to the quantum coherent effects studied recently 
in the photosynthetic complexes Ba¬ 
it is important to note, that in many systems the occurrence 
of the ST requires not only the overlapping of the neighbor¬ 
ing resonances, but also a delocalization of the eigenfunctions 
which are involved in the ST. This is especially important for 
those systems in which the initial wave function does not suffi¬ 
ciently overlap with those eigenfunctions which are responsible 
for the SR (see below). 

The ST has been well-studied in random matrix theory 
BEiijaiiiaiigin], nuclear physics [1811191120] . microwave bil¬ 
liards om m], and in paradigmatic models of coherent quan¬ 
tum transport [231 M [25]. In particular, it was shown in 
[ 231123 ] that, for a realistic model for quantum transport, max¬ 
imum transmission is achieved at the ST. In biological systems, 
the ST has previously been discussed in [TTl fT2~t fSUt [ 27 : . (See 
also references therein.) 

In this paper, we numerically study the possible role of the 
ST in the NPQ mechanism in the LHCs of photosystem II 
(PSII). For concreteness, we apply our approach to the NPQ 
protocol discussed in [28]. According to [28], the NPQ mech- 
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anism is associated with charge transfer quenching. It takes 
place in part in the LHC CP29 minor LHC, which is located 
(with other minor LHCs) between the RC and the peripheral 
LHCII complex [25]. In the NPQ mechanism, discussed in [281 . 
the following basic components are involved: a dimer based on 
two chlorophylls with their excited states, Chla% and Chlbt,- 
the heterodimer, ( Chla 5 — Zea)*, where the carotenoid zeax- 
anthin, Zea, is formed in the NPQ regime; the charge transfer 
state (CTS) of this heterodimer, (Chla§ — Zea + )*; and the 
damaging channel, which is responsible for production of the 
singlet oxygen and other damaging products. Our model (see 
also 0) includes five discrete energy levels (sites), each for 
a mentioned above component. Each of two discrete levels, 
associated with the CTS and with the damaging channel, is 
connected to its corresponding independent sink, which is de¬ 
scribed by the continuum electron spectrum. Both these sinks 
are responsible for the quasi-reversible chemical reactions. Ac¬ 
cording to [25], in the NPQ regime, the extra sunlight energy, 
that is absorbed by LHC, is transfered to the CTS, and then 
it partly dissipates into the protein environment. 

In our model, all five discrete states and two sinks are em¬ 
bedded in the protein environment, which we model by a ran¬ 
dom telegraph process (RTP). In [2], our model was used to 
simulate numerically the efficiency of the NPQ regime, at vari¬ 
ous parameters of the system. In the present paper, we extend 
the application of our model to describe the possible role of the 
ST in the NPQ regime. Our mathematical approach is based 
on a non-Hermitian Hamiltonian which naturally occurs and 
is a standard mathematical tool used for studying the ST in 
many physical and biological systems. We demonstrate that, 
for reasonable conditions, the NPQ mechanism can be asso¬ 
ciated with the ST in the charge transfer channel. We also 
show that the ST can optimize the performance of the NPQ 
mechanism. 


Mathematical formulation of the model 

Our model includes five discrete excitonic energy states and 
two sinks (see Fig. ©. Below, we use the notation and the 
basic NPQ components similar to those used in [2] [25]. The 
values of the parameters are varied in order to demonstrate 
the different dynamical regimes of the NPQ mechanism and 
their relation to the ST. We use the following notation. The 
discrete electron state, |1) = \Chla% ), is the excited electron 
state of Chla 5 . The discrete electron state, |2) = \Chlb%), is 
the excited electron state of Chlb$. So, in our model a dimer 
based on two chlorophylls, Chlas and Chlb$, each in an ex¬ 
cited electron state, participates in the NPQ process, as in [25] . 
The discrete state, |3) = \(Chla$ — Zea)*), is the heterodimer 
excited state. The discrete state, |4) = |(a~ — Zea + )*), is 
the CTS of the heterodimer. The state, |0) = \damage), is a 
discrete part of the damaging channel. The sink, |Sb), is the 
continuum part of the damaging channel. The sink, IS 4 ), is 
the continuum part of the CTS (channel), |4). All matrix ele¬ 
ments, Vmn, in the Hamiltonian © describe the interactions 
between the discrete excitonic states. Each sink is character¬ 
ized by two parameters: the rate, To ^ 4 ), of ET into the sink, 
|Sb) (IS4)), from the corresponding attached discrete state, | 0 ) 
(|4)), and the efficiency (cumulative time-dependent probabil¬ 
ity), r]o(t) ( r/ 4 (t )), for the exciton to be absorbed by the cor¬ 
responding sink. Note that both Po and T 4 characterize only 
the rates of destruction of the exciton in the CP29 LHC, and 
they do not describe any subsequent chemical reactions that 
take place in the damaging and the charge transfer channels. 
In this sense, our model describes only the primary NPQ pro¬ 
cesses in the ET, and it does not describe the processes which 
occur in both sinks. The latter occur on relatively large time¬ 


scales, and require a detailed knowledge of the structures of 
the sinks, and additional methods for their analysis. 

According to [2], the quantum dynamics of the ET can be 
described by the following effective non-Hermitian Hamilto¬ 
nian, T-L = H — zW, where, 

U = y~[gn|n)(n| + y~] Vmn\m) (n\, m,n = 0,1,...,4, [1] 

n m^n 

is the dressed Hamiltonian, and 

t 2 ] 

n= 0,4 

In ©, £ n is the renormalized energy of the state, \n), and 
the parameter, T„, in © is the tunneling rate to the n-th sink. 
Below, only To and T 4 are taken into account. 

The dynamics of the system can be described by the 
Liouville-von Neumann equation, 

p = i[p,H]-{W,p}, [3] 

where {W, p } = Wp + pW. 

We define the ET efficiency of tunneling to all N sinks as, 

z?(t) = 1 - Tr(p(f)) =* f Tr{W,p(r)}dr. [4] 

Jo 

This can be expressed as the sum of time-integrated probabili¬ 
ties of trapping an electron into the n-th sink, r/{t) = r] n (t) 
[30] [3T]. Here, 

r/n(t) =T n / pnn(T)dr, [5] 

Jo 

is the efficiency of the n-th sink. In particular, for our model, 
presented in Fig. [I] complete suppression of the damaging 
channel at all times occurs if r]o(t) = 0 . 

The eigenvalue problem for the non-Hermitian Hamiltonian. 

The bi-orthogonal eigenstates of the effective non-Hermitian 
Hamiltonian, H, provide a convenient basis in which the eigen¬ 
value problem can be formulated and resolved [T], [TD]- (See also 
references therein.) Let \>pa) be the right eigenstates of H with 
complex eigenvalues, E a , 

H\ip a ) = E a \*p a ), a = 1,2... 5. [ 6 ] 

Further it is convenient to set E a = £ a — iT a , where, E a and 
T a (decay width), are the real and the negative imaginary 
parts of the complex eigenvalue. 

By manipulating the tunneling rates into the sinks, one can 
modify the dynamics of the whole system. As will be demon¬ 
strated below, under some conditions, quantum coherent es¬ 
cape of the exciton in the charge transfer channel occurs. This 
escape is associated with the ST in the sink, IS 4 ). A similar 
ST can be realized, under some conditions, in the damaging 
channel (sink |S'o)). 

Noise-assisted ET in the NPQ regime. In the presence of noisy 
environment, the evolution of the system can be described by 
the following effective non-Hermitian Hamiltonian [2], 

n to t = n-iw + v(t), [7] 

where, V(t) = A m „(t)|m)(n|, m,n = 0 , 1 ,... 4, 

A mn{t), describes the influence of the protein noisy environ¬ 
ment. The diagonal matrix elements of noise, A nn (t), are re¬ 
sponsible for decoherence. The off-diagonal matrix elements, 
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A mn{t), lead to direct relaxation processes. In what follows, 
we restrict ourselves to the diagonal noise effects. (See also 
|32l 1331 134: . and references therein.) Then, one can write, 
A mn{t) = ^nSmn^{t), where A n is the coupling constant at 
site, n, and /(f) is a random process. 

We describe the protein noisy environment by an RTP, /(f), 
with the following properties, 

m) = o, is] 

= [ 9 ] 

where, y(f — t') = a 2 e~ 2l ^ t ~ t \ is the correlation function of 
noise. 

The evolution of the average components of the density ma¬ 
trix is described by the following system of ordinary differen¬ 
tial equations [2]: 

j t {p)=- (w, {p)\ - [io] 

^<P 5 > = i[(p ( ),H] - {W, (/> - iB(p) - 2 7 (pt)[ll] 

where (p 5 ) = (/p)/o-, B = J2m,n( d ™ ~ d n )\m){n\, and we set 
dm = A m cr. The average, (...), is taken over the random pro¬ 
cess. In deriving Eqs. ffiB and CD, the approach developed 
in GSHMIII] for the RTP was used. 

Employing Eqs. tna and HI one can show that the follow¬ 
ing normalization condition is satisfied, 

4 

^(Pnn(t))+Po(t) + P4(t) = l. [ 12 1 

71 = 0 

Eq. m requires that the total probability of finding the ex- 
citon among the five discrete levels and in two sinks is unity 
for all times. 


Results of numerical simulations 

Here we examine numerically the possible role of the ST in 
the NPQ regime. The results of our numerical simulations 
are obtained by solving Eqs. m and CD for the average 
components of the density matrix. Below, in numerical sim¬ 
ulations, we choose h = 1. Then, the values of parameters in 
energy units are measured in ps _1 (lps -1 « 0.66meV). Time 
is measured in ps. 

To simplify the presentation of our results, some charac¬ 
teristic parameters were fixed. The following values of these 
parameters were chosen: The energy of the excited state of 
Chla%, ei = GOps^ 1 ~ 39.6meV; The energy of the excited 
state of Chlb%, £2 = 210ps _1 « 138.6meV; The energy of 
the excited state of the heterodimer, (Chla 5 — Zea)*, £3 = 
45ps -1 « 29.7meV; The energy of the CTS of the heterodimer, 
(ChlaJ — Zea + )*, £4 = 30ps _1 « 19.8meV. The following 
values of the matrix elements of the dipole-dipole interactions 
between the sites were fixed: V 10 = V 20 = 30ps _1 ~ 19.8meV, 
V 12 = 15ps _1 « 9.9meV, V 13 = V 34 = 25ps _1 « 16.5meV. 
The values of the coupling constants were chosen: d\ = £ 1 , 
d 2 = £ 2 , do = £ 3 , dr = £4. These values of the coupling 
constants provide the strongest (resonant) noise on the cor¬ 
responding sites (see below). The correlation decay rate was 
hxed: 7 = lOps^ 1 . The time duration for simulation of the 
ET efficiencies, po(r) and 74 (r), was fixed at r = 5ps. The 
parameters: £0 (characterizing the energy level of the dam¬ 
aging state, |0}); the coupling constant, do, of the state, |0), 
with noise; and the ET rates, Po and P 4 , were varied in some 


simulations. The initially populated state was chosen to be 
the excited state of Chlat,. 

In Figs. [2ji and Hb we computed the efficiencies of both 
sinks, 770 , 4 ( 1 "), in the absence of noise. The energy of the dam¬ 
aging state, | 0 ), was chosen: £0 = 0 , close enough to initially 
populated level, £1, of Chla 5. As one can see from Fig. Ht, 
when only one sink, IS 4 ) (related to the CTS) is open (green 
curve, Fo = 0), the efficiency of the NPQ process experiences 
a maximum [rf^ ax « 97%), at the ET rate, F° « 72ps _1 . The 
important observation is that when F 4 increases, the efficiency, 
774 , decreases. As will be discussed below, this dependence of 
the efficiency, 74 on T 4 is caused by the ST. When two sinks 
are open (Fo = 2ps _1 ), the efficiency of the NPQ process 
(blue curve in Fig. Hr) still has a maximum as a function of 
the ET rate, T 4 . In this case, rf± a:L « 82% and F° a: 72ps _1 . 
This decrease of the NPQ efficiency is caused by the damaging 
channel, ISb), during r = 5ps (red curve in Fig. Hi). I 11 Fig. 
H, the opposite case is considered, when the ET rate to the 
CTS (sink IS 4 )) is fixed (F 4 = lOps -1 ), and both efficiencies, 
770 and 74 , are simulated as functions of the ET rate, To, to 
the damaging sink. Similar to the results presented in Fig. 
H, in Fig. H one can observe the ST in the damaging sink, 
|S'o). In this case, the maximal efficiency is relatively small, 
i )™ 1 « 70%, because the ET rate to the CTS is relatively 
high (F 4 = 10 ps _1 ). 

In Figs. H> and Hb the results of numerical simulations 
are presented for the NPQ efficiency, 774 ^ 4 ) and for the ef¬ 
ficiency of the damaging channel, 770 ^ 4 ), for hxed ET rate, 
To = 2ps _1 . So, two sinks are open in this case. The two main 
differences from the previous results presented in Figs. H 
and Et>, are caused by the following two modifications. First, 
in Fig. [ 3 ] 3 , the energy of the damaging state was chosen, 
£0 = —90ps _1 « — 59.4meV, which is significantly below the 
value, £ = 0, which was used in Figs. [5Jt and[5j7 (and in Fig. 
Hi). Second, the influence of noise, produced by the protein 
environment, is included. 

In our model, noise is generated by the RTP, and it has 
two characteristic parameters: the amplitude of the correla¬ 
tion function, a 2 , and the correlation time, r c = 2 / 7 . Also, 
the influence of noise on the CP29 LHC is characterized by the 
interaction constants (matrix elements of “site-protein” inter¬ 
actions), A n (n = 0 , 1 ,.., 4). Then, instead of parameters, A n , 
the renormalized parameters, d„ = \ n cr, occur in our model. 
The important property of the influence of noise on the chloro¬ 
phyll dimer, discussed in [3H][35], is the following. Under the 
“resonant conditions”, £„ — £ m = d n — d m , the influence of 
noise on the dimer, n — m, becomes maximal in the sense that 
the ET rate has a resonant peak. 

In Fig. Hi, the damaging level (£0 = 0) is chosen close to the 
energy level of Chlat,. The efficiencies of sinks, 774 and 770 , are 
indicated by the blue and red curves, respectively. Blue and 
red solid curves correspond to the resonant noise (do = 0 ) on 
the dimer, 1 — 0. Blue and red dashed and blue and red dash- 
dot curves correspond to non-resonant noise on the dimer, 
1 — 0. Green and yellow curves correspond to the absence of 
noise. As one can see, the maximum of the NPQ efficiency 
is, ri™ ax « 85%, almost independent of noise, in the vicinity 
of r!j « 72ps _1 . This behavior of 774 , in the vicinity of F°, 
can be explained as follows. Because the damaging sink is 
close (in energy) to the Chlat, energy level, the damaging sink 
accumulates during r = 5ps the amount of probability which 
practically is independent of the amplitude of noise. 

Different dependence of the NPQ efficfency on nofse, in the 
vicinity of its maximum, occurs when the damaging energy 
level, £0 = —90ps _1 « 59.4meV, is relatively far from the en¬ 
ergy level of the Chla%. As one can see from Fig. Hb the res¬ 
onant noise (do = £ 0 ), acting on the dimer, 1 — 0 , significantly 
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degrades the NPQ performance in the vicinity of ri™ ax « 88% 
(blue dashed curve). The reason is that the resonant noise 
acting on the 1 — 0 transition, helps to accumulate extra prob¬ 
ability in the damaging sink. When noise is non-resonant on 
the 1 — 0 transition (blue curve) the NPQ efficiency improves 
up to, T]™ ax m 95%, and it approaches the NPQ efficiency 
when noise is absent (green curve). 

The superradiance transition in relation to the l\IPQ. When 
noise is absent, the ST, in relation to the NPQ mechanism, is 
easy to explain in terms of the effective non-Hermitian Hamil¬ 
tonian, H = Tl — iW, where the Hermitian operator, H , is 
defined by Eq. ©, and the Hermitian operator, W, is defined 
by Eq. In Eq. @, only Po and P 4 differ from zero. The 
eigenvalue problem for Iri is formulated in Eq. ©. In Fig. 
EH, the blue curve corresponds to the negative imaginary part 
(decay width) of the eigenvalue, Y a , with T max, as a function 
of T 4 (blue curve). The red curve corresponds to the average 
decay width, T a „, of 4 eigenstates with smallest values of T a , 
as a function of T 4 . The energy level of the damaging state 
was chosen, so = — 90ps _1 ~ —59.4meV. As one can see, in 
the vicinity of PJ « 84ps _1 , the red curve experiences a maxi¬ 
mum. At this point, a segregation of the decay widths begins. 
Namely, when T 4 grows, the red and blue curves start to sep¬ 
arate. So, the decay width of one (superradiant) eigenstate 
continue to grow with P 4 , and the decay widths of all other 
(subradient) eigenstates decreases. This segregation behavior 
of the decay widths is one of the characteristic properties of 
the ST [3 [TO]. As one can see from Fig. EH (green curve), 
in the vicinity of the ST (rj « 84PS” 1 ), the NPQ efficiency 
experiences its maximum behavior, with r/^ lax « 95%. 

Similar ST behavior occurs in Fig. EH, when T 4 = 0, and 
the dependence of widths on tunneling rate into the sink |5o) 
is presented. The insert shows the detailed behavior of the red 
curve. In this case, the ST occurs in the damaging sink, |S'o). 

Overlapping of resonances. In Figs. EH and EH, we show the 
dependences on T 4 of the decay widths, T a , and the real 
parts of the eigenvalues, £ a . Different colors correspond to 
five different eigenstates. As one can see, at the superradi¬ 
ance transition (rj « 84ps _1 ps), the superradiant state, with 
largest decay width (blue curve), begins to overlap (interact) 
with its neighboring eigenstate (red curve). Indeed, at this 
point the spacing between these two states (resonances) is 
approximately equal to the half-sum of their decay widths: 
Sbiue — £red « T blue — T rec j. Then, these resonances cannot be 
considered to be isolated, and they both coherently contribute 
to the ET to the corresponding sink (continuum) [3 [10]. In 
our case, only two resonances overlap at TJ. When P 4 grows, 
the decay width of the superradiant state grows, and finally 
the superradiant resonance overlaps with all other resonances. 
(See Figs. EH andEH-) But this does not mean that the super¬ 
radiant state will continue to contribute (for large ET rates, 
P 4 ) to the ET into the continuum (sink). The localization 
properties of the eigenstates become important. (See the next 
sub-section.) Note also, that usually many interacting reso¬ 
nances may overlap and contribute to the ST El ED]. 

The trajectories of the eigenenergies. When discussing the 
ST, the diagram of the trajectories of the eigenenergies in the 
complex plane (3 = Sr(-E), £ = 3t(£/)) can produce additional 
insight to the features of the ST mm- We present this di¬ 
agram in Fig. [6] for five eigenenergies (indicated by different 
colors) parametrized by the ET rate, T 4 £ [0,10 4 ], and for the 
fixed value of To = 2. The superradiant eigenenergy is indi¬ 
cated by the blue curve. Its decay width increases when T 4 


increases. The eigenenergies of four other eigenstates become 
trapped, as T 4 increases. Their decay widths are bounded, 
and experience a maximum as a function of IV The super¬ 
radiant eigenstate (blue curve) overlaps with the neighboring 
eigenstate, indicated by the red color, at TJ « 84ps _1 . All 
other neighboring resonances do not overlap. The detailed 
behaviors of the trajectories of two eigenenergies, located in 
the upper left and the upper right, are shown in the insets, 
(a) and (b). The insert (c) demonstrates the behavior of the 
eigenenergies for large values of T 4 . The following parameters 
were chosen: so = — 90ps _1 « — 59.4meV and To = 2ps~ x . 

Localization of the eigenfunctions. To study numerically the 
localization properties of the eigenstates, we use the partici¬ 
pation ratio (PR) [10] of the eigenfunction, \ip a ), defined as, 

PRa = ( ^ |<n|V’a>| 4 ) _1 . [13] 

n 

Its value varies from 1, for a fully localized eigenstate, to 5, 
for a fully delocalized eigenstate. In Fig. 0 the PRs for all 
five eigenfunctions are presented, as functions of 74 . As one 
can see, the superradiance state (blue curve) is partly delocal¬ 
ized for relatively small values of T 4 < TJ « 84ps _1 , and it 
becomes localized at the site, |4), for large values of T 4 . The 
sub-radiant eigenstates also experience reorganization in the 
vicinity of 1% < TJ « 84ps _1 , but they remain partly delo¬ 
calized for large values of T 4 . These sub-radiant eigenstates 
contribute to the NPQ mechanism for large values of 1%. 

In Fig. EH, we demonstrate the localization properties of 
eigenstates, for a moderate value of To = 2 ps _1 , and for a 
relatively large value of T 4 = 10 3 ps _1 . The colors correspond 
to those used in Figs. [6]and[3 On the vertical axis, the prob¬ 
ability, P a (n ) = \(n\tl> a }\ 2 , is presented, of the population of 
the site for a given eigenfunction, |i p a ). Different colors indi¬ 
cate different eigenfunctions. As one can see, the superradiant 
eigenstate, indicated by the blue curve, is localized mainly on 
the CTS, |4). This is the main reason why this eigenstate, 
which has the largest decay width, does not contribute signifi¬ 
cantly to the ST, or to the NPQ dynamics. (See green curve in 
Fig. EH.) Indeed, the initial condition corresponds to the pop¬ 
ulation of the Chla% (site |1)). Then, at large values of the 
ET rate, T 4 , this initially populated state, 11), only weakly 
overlaps with the superradiant eigenstate. So, during the fi¬ 
nite time, r = 5ps (as in Figs. EH andEH)), the efficiency of 
the NPQ, 74 (r = 5ps) is relatively small. The NPQ efficiency 
remains relatively small, for large values of T 4 , in the presence 
of noise. (See, Figs. EH andEH-) 

In Fig. EH, we show the localization properties of the 
superradiant eigenstate (with the largest decay width), for 
r 2 = 2ps _1 , and for four values of the ET rates, IT. As 
one can see, for relatively the small value, T 4 = 10 ps _1 (green 
color), the superradiant eigenstate is inhomogeneously delo¬ 
calized over all five sites. The value T 4 = lOOps ^ 1 (blue color) 
corresponds to the vicinity of the ST, and to the best perfor¬ 
mance of the NPQ mechanism. The superradiant eigenstate is 
still inhomogeneously delocalized over all five sites. When T 4 
increases (T 4 = 200 ps _1 , orange color), the localization of the 
superradiant eigenstate on the CTS, |4), increases. For large 
values of T 4 = 10 3 ps ^ 1 (red color), the superradiant eigenstate 
is mainly localized on the CTS (site |4)). 

Conclusion 

Our main result is the demonstration of the relation be¬ 
tween the two well-known phenomena: the nonphotochemical 
quenching (NPQ) mechanism in photosynthetic organisms and 
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the superradiance transition (ST), that occurs in many physi¬ 
cal and biological systems. For concreteness, we analyzed this 
relation by using a relatively simple model for the NPQ by 
the charge transfer state (CTS), in the CP29 minor LHC of 
the PSII. The NPQ protocol we used, is mainly based on the 
experimental results discussed in 55] , 

In the absence of noise, our model is described by an ef¬ 
fective non-Hermitian Hamiltonian, which is represented by a 
non-Hermitian 5x5 matrix. The eigenenergies of this matrix 
are complex numbers. According to the accepted terminol¬ 
ogy, their real parts represent the positions of the resonances 
(poles), in the complex energy plane. The absolute values of 
their imaginary parts are the decay widths of the resonances. 
The ST occurs when these resonances begin to interact (over¬ 
lap). In this case, the superradiant eigenstates accumulate 
large decay widths. So, these superradiant states contribute 
coherently to the radiation into the continuum spectra (sinks). 
All other eigenstates (sub-radiant) have small decay widths, 
and are long-lived. However, for some range of parameters, the 
superradiant eigenstates become localized at the sites which do 
not overlap with those sites which are initially populated. In 
this case, the superradiant eigenstates practically do not con¬ 
tribute to the ST (in spite of the fact that these eigenstates 
have the largest decay widths). This competition between the 
decay widths and the delocalization of the superradiant eigen¬ 
states is well-know in the field of the SR. Since our model 
includes two sinks, two STs occur for different values of pa¬ 
rameters. One ST occurs in the damaging channel, and the 
other occurs in the charge transfer channel. This latter one 
contributes to the NPQ mechanism. 

The efficiency of the NPQ by the CTS mainly depends on 
two factors: the ET rate to the sink connected to the site 
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of the CTS, and on the probability accumulated at the ex- 
citonic site, where the CTS is localized. When the ET rate 
to the sink associated with the CTS is small, the NPQ effi¬ 
ciency is also small. However, when this ET rate is very large, 
the superradiant eigenstate becomes strongly localized at the 
excitonic CTS, and it practically does not overlap with the ini¬ 
tially populated chlorophyll dimer. Then, the NPQ efficiency 
is also suppressed in this case. Only for the intermediate val¬ 
ues of the ET rate in the sink associated with the CTS, does 
the NPQ efficiency approaches its maximum. 

We demonstrated that the maximum of the NPQ efficiency 
is caused by the ST to the sink associated with the CTS. 
However, experimental verifications are required in order to 
determine whether or not the NPQ regime is associated with 
the ST transition for real photosynthetic complexes. Indeed, 
it can happen that, in the photosynthetic apparatus, the NPQ 
regime occurs in the “non-optimal” region of parameters, and 
it could be independent of the ST. 

The results obtained in this paper can easily be generalized 
to more complicated networks of interacting chlorophyll and 
carotenoid light-sensitive molecules, and for more complicated 
random protein environments. The experimental verification 
of the relation between the NPQ and the ST phenomena would 
be of significant interest for both better understanding of the 
NPQ mechanisms and for bio- and nano-engineering. 
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Fig. 1. Schematic of our model consisting of five discrete states, |n), (n = 0,_, 4), 

and two independent sink reservoirs, \So) (connected to the damaging state), and IS 4 ) (con¬ 
nected to the charge transfer state). The ellipse indicates the dimer based on excited states of 
two chlorophylls. The blue dashed lines indicate non-zero matrix elements used in numerical 
simulations. 



Fig. 2. (Color online) The efficiencies of two sinks in the absence of noise ( £0 — 0)- ( a ) The efficiencies of sinks, 774 (blue and green curves) and 770 (red curve) vs 
^r 4 , for To =0 (green curve) and for Tq = 2 (blue and red curves); (b) The efficiencies of sinks, 774 (blue curve) and 770 (red curve) vs. ^/Fq, for F 4 = 10. 
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Fig. 3. (Color online) The efficiencies of two sinks in the presence of noise: 774 (blue curve) and 770 (red curve) vs. y/T 4. (a) £0 = 0: solid curves: resonance noise for 
dimer 1 — 0 (do = 0); dashed curves: do = —30; dash-dot curves: do = 30. (b) £0 — —90: solid curves (do = 0); resonance noise, do = —90 (dashed curves). 
r 0 = 2 . Green and orange curves demonstrate the efficiencies of sinks in the absence of noise. 




Fig. 5. (Color online) (a) Decay widths as functions of the rate F4. (b) 8 = ?R.(E) VSF4. The blue curves correspond to the largest decay width. (£q = — 90, Tq = 2 ). 
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Fig. 6. (Color online) The trajectories of five eigenvalues in the complex eigenenergy plane (T = ^s{E), £ = $?(£/)), as F 4 increases from 0 to 10000. The superradiant 
eigenstate (blue curve) is "unstable” (its decay width, |T|, increases when r .4 increases). All other four eigenstates become trapped. (Their decay widths are bounded above.) 
The insets, (a) and (b), show the detailed behavior of two eigenvalues located in the upper left and the upper right. The insert (c) demonstrates the behavior of the eigenvalues 
for large values of (£q — — 90, Fq = 2 ). 



Fig. 7. (Color online) The participation ratios of the eigenstates of the effective non-Hermitian Hamiltonian, l~i, as functions of (£q = —90, Tq = 2). 
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Fig. 8. (Color online) (a) Localization of the eigenstates (r .4 = 10 3 ). (b) Localization of the eigenstate with the largest decay width: T 4 = 10 (green), T 4 = 100 
(blue), r 4 = 200 (orange), T 4 = 10 3 (red). (£q = — 90, Tq = 2). 
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